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SOFT X-RAY EMISSION LINES IN THE EARLY AFTERGLOW OF GAMMA-RAY BURSTS 
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ABSTRACT 

We compute the luminosity of Ka emission lines produced by astrophysically abundant elements in the soft 
X-ray spectra of the early afterglow of gamma-ray bursts. We find that the detection of these lines can be a 
diagnostic for the geometrical set-up of the reprocessing material. In particular we can distinguish between a 
"geometry dominated" model, in which the line emission is coming from an extended region and its duration 
arises from light-travel-time effects and an "engine dominated" model, where the line emitting gas is in a smaller 
region, irradiated for a longer period. These lines therefore offer clues to the dynamics and time-scale of the 
explosion leading to a gamma-ray burst. 

Subject headings: gamma rays: bursts — line: formation — radiation mechanisms: non-thermal 



L INTRODUCTION 

Narrow emission features superimposed on the early X-ray 
afterglow of gamma-ray bursts (GRBs) seem to be common 
[Piro et al. 1999, AntonelH et al. 2000, Piro et al. 2000 (here- 
after POO), Reeves et al. 2002 (hereafter R02)], even though 
a single strong evidence of their significance is still lacking. 
These features are observed starting from several hours after the 
burst (due to instrumental limitations) up to a couple of days. 
Some have signs of variability (e.g. R02), others are constant 
for the total duration of the observation (Antonelli et al. 2000; 
POO). Understanding how and where they are produced is fun- 
damental for identifying the nature of the burst progenitor In 
fact the line, being observed at its rest wavelength, has to be 
produced through reprocessing of the burst radiation, and car- 
ries therefore information about the geometry and structure of 
the material surrounding the burst (Lazzati, Campana & Ghis- 
ellini 1999). It is commonly accepted that the presence of metal 
lines itself strongly favors hypernova models (Woosley 1993; 
Paczynski 1998) and can be used to rule out compact merg- 
ers (Eichler et al. 1989) as burst progenitors. Nevertheless, 
whether the lines can be produced in a standard hypernova sce- 
nario or require an explosive event or supernova that occurred 
prior to the GRB (see e.g. Vietri & Stella 1998) is still a matter 
of open debate. 

We define as Engine Dominated (ED) the models in which 
the lines are produced in a standard - single step - hypernova 
explosion. In these models, the line is created by reprocessing 
from material very close to the explosion site (R ~ 10^^ cm). 
The ionizing continuum in this case can not be attributed to the 
burst or afterglow itself, since the photons are radiated at larger 
distances. The ionizing continuum is then believed to be pro- 
vided by a long lasting engine (Rees & Meszaros 2000) or by 
magnetic energy stored in a plasma bubble emerging after the 
jet has crossed the stellar progenitor (Meszaros & Rees 2001). 
In ED models the duration of line emission is related to the du- 
ration of the the ionizing continuum. 

We define as Geometry Dominated (GD) the models in which 
the reprocessing material is located at a large enough distance, 
R, to be illuminated by the burst and afterglow photons. In 
these models the duration of the line emission is set by the size 
of the reprocessor This reprocessing material has to be com- 



pact and metal enriched, similar to a supernova remnant (SNR), 
as naturally predicted in the supranova model (Vietri & Stella 
1998). In both ED and GD models the line is supposed to be 
produced by reflection off a slab of optically thick material (see 
also Bottcher 2000 and Bottcher & Fryer 2001 for an alternative 
GD model). 

In this letter we compute the luminosity of Ka lines from the 
elements O, Ne, Mg, Si, S, Ar and Ca (which we will refer to 
as "light elements" hereafter), and compare them to the iron 
one. We find that the detection of soft X-ray (~ [0.5-5] keV) 
features, possibly made in the afterglow of GRB 011211 (R02) 
strongly favors GD scenarios for the production of metal lines 
in GRB afterglows. The letter is organized as follows. In §2 
we compute the line luminosities in reflection from an a priori 
point of view. We compute the line flux in GD and ED sce- 
nario in §3 and we discuss the comparison of predictions with 
observational constraints in §4. 

2. LINE LUMINOSITIES 

The luminosity of the lines has been computed according to 
the equation: 

Lx = ^ (1) 

where the subscript x indicates the particular ion, «x is the ion 
density, ex is the Ka emission line frequency ex = hvx of the 
H-like ion, S is the emitting surface, A is the depth of the emit- 
ting layer and fen, is the time scale of emission of a photon by a 
single ion. 

The time scale of photon emission can be set either by the 
ionization time scale corrected for the Auger effect (if the ion- 
ization parameter is low) or by the recombination time scale, for 
highly ionized plasma. Here, we compute the emission time as 
max(f,-ec, fionMwjx.z}), where Am{x,z} is the photoelectric yield 
for the element X with charge Z (Kaastra & Mewe 1993). The 
ionization time scale fion is computed as (Lazzati et al. 2001b): 
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where L(j/) = Lqv "is the ionizing spectrum (note that Lq does 
not have the dimensions of a luminosity), ax the threshold cross 
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section of the element X (atomic data are taken from Verner 
& Yakovlev 1995) and R the distance of the reprocessing ma- 
terial from the ionizing flux source. The recombination time 
scale free has been computed with the analytic approximations 
and tables of Verner & Ferland (1996). A constant temperature 
T = 4x 10^ K has been assumed. Both the ionization and re- 
combination time scales are computed for the H-like ion. The 
charge of the atom is then computed as' Zx = (fion/fi-ec)^''^- 

The depth of the emitting layer A is assumed to be the depth 
that the ionizing photons can reach and/or a line photon can es- 
cape without a shift in frequency. We therefore compute it as 
the minimum between the plasma opacity at resonance for the 
ion X and the maximum Thompson depth a photon can travel 
before being shifted by Thomson scattering. We therefore have 
A = min(A7-, Ax), where: 



Y 



-\ -1 



f{ion,Y} +?{rec,Y} 



and 



At — 



tie CTt 



(3) 



(4) 



where the sum in Eq. ^ is extended over the astrophysically 
relevant elements H, He, C, N, O, Ne, Mg, Si, S, Ar, Ca, Fe (el- 
emental abundances taken from Anders & Grevesse 1989), and 
o'f(j^x) is the cross section of element Y computed at the thresh- 
old frequency of element X. In the numerator of Eq. ^ N^dT) 
is the maximum number of scatterings a photon can make be- 
fore the line is smeared or shifted. If e ~ kT, line broadening is 
more important than shifting, and we have a^/e ~ ^JkT /{nie c^) 
for a single scattering. Requiring a limit of (Xj/e < 1/3 in or- 
der for the line to be detectable as a narrow feature, we have 
Tmax = V^Z ~ [nieC^ /i9kT)]^l'^ « 4. Note however that, for 
the considered abundant elements, Eq. ^ is relevant only for 
very high ionization parameters (i.e. in the decreasing tails of 
Fig. |l]) and therefore its uncertainty is not necessarily reflected 
in the maximum efficiency of line reprocessing. 

In Fig. |l] and ^ (lower panel) we show the results of our cal- 
culations for O, Ne, Mg, Si, S, Ar, Ca and Fe as a function of the 
ionization parameter ^ = 47rF[i_i()]/n, where F[i-io] is the ioniz- 
ing flux in the energy range [1-10] keV and n the particle den- 
sity. The figures show that, in optimal conditions, ^ 1 % of the 
ionizing flux is reprocessed in each line. We also show that in- 
creasing the metallicity does not increase this fraction, but helps 
keep it constant for high ionization parameters. The accuracy 
of this line strength estimates have been checked against the re- 
sults of numerical simulations for Mg, Si and Fe (Ross 1979; 
Ross & Fabian 1983; BaUantyne, Fabian & Ross 2002). We 
find that our results are in agreement within a factor ~ 2 with 
line luminosities obtained by fitting Gaussian profiles to the nu- 
merical spectra. The largest deviations are observed close to the 
minimum value of the efficiency (in the strong Auger regime), 
where our semi-analytical approach yields larger luminosities 
with respect to the full numerical treatment. This comparison 
has been performed with Mg, Si and Fe lines, since the numer- 
ical code does not include S, Ar and Ca. 



In order to apply the results of the previous section to various 
GRB models, we need to transform the ratios shown in Fig. [l| 
and ^ into observed fluxes and/or luminosities. Even though 
in spectroscopy it is customary to compute the line equivalent 
width, in the case of GRBs it is assumed that the observed con- 
tinuum is not related to the ionizing continuum responsible for 
producing the lines (which can be the afterglow itself at earlier 
times in GD models or a completely different component in ED 
models), and so it is more useful to directly compare the line 
fluxes (Lazzati et al. 1999). 

3.1. Geometry dominated models 

In the case of GD models, the observed line flux has to be 
corrected by a geometric factor This factor takes into ac- 
count that while the emitting material is illuminated for a time 
fill, the line is observed, at infinity, for a time fobs ^ Rsnr/c- A 
precise computation of ^ = fm/fobs is complex, since the re- 
processing material is illuminated by the burst and early af- 
terglow radiation, and it is likely that the properties of the 
fireball towards the reprocessing material are different from 
those along the line of sight (Rossi, Lazzati & Rees 2002). 
If one assumes that the fireball decelerates inside the remnant 
and that the density of the external medium in the immedi- 
ate vicinity of the progenitor scales with the remnant radius^ 
as n = 10'^n8(^SNR/10'^)"-' cm"^, the time the firebafl takes to 
overtake the remnant is proportional to the radius of the rem- 
nant itself, and: 
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where fcRB is the burst duration, and £51 is the isotropic equiv- 
alent energy output of the burst in units of 10^' erg in the di- 
rection of the reprocessing material. We have used Eq. 15 of 
Panaitescu & Kumar (2000) to evaluate the numerical factor. 

On the other hand, in GD models the ionizing continuum 
is not directly observed, and can then be assumed to be much 
larger than the afterglow luminosity at the time the line is de- 
tected. The ratio 77cont of reflected to incident continuum de- 
pends on the ionization parameter, ranging from < 0.1 for 
^ < 10^ to ~ 1 for ^ > 10^ (see the upper panel of Fig. For 
reasonable ionization parameters, one can then assume an ion- 
izing luminosity Lion < i'After/(C'7cont) '^-^ 2 X 10'*^ erg s"', where 
After is the afterglow [1-10] keV luminosity at the time at which 
the lines are observed. The apparent contradiction of having an 
ionizing continuum brighter than the afterglow but not observed 
in the afterglow itself, is solved if we consider the geometry of 
the system (see e.g. panels a and b of Fig. 2 in Vietri et al. 
2001). In fact, the fireball emission towards the reprocessing 
material is not necessarily related to the one in the direction of 
the observer Moreover, line photons produced at time f are ob- 
served at a later time due to the longer travel thay have to make. 
At this later time the afterglow emission has decreased. 

Finally, we must take into account that while the ionizing 
continuum may be geometrically beamed, the line emission is 
nearly isotropic, so that the observed flux must be corrected 
by a factor n/A-n ^ 0.25, where the numerical factor has been 
computed for a jet opening angle of 6* = 40°. This open- 
ing angle is one order of magnitude larger than the opening 

'The exponent 2/5 has been derived by interpolation of the numerical recombination rates of Shull & van Steenberg (1982). 

^Note that the density is not supposed to scale with the radius. We consider a uniform density inside the remnant, scaling it with the remnant radius in order to 
have a constant total mass inside it. 
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angle measured in GRB 991216 (Frail et al. 2001). Since 
J7 /47r = 1 - cos 0y oc 6* J , a ten times smaller opening angle would 
make the lines two orders of magnitude fainter and thus unde- 
tectable. We here assume this larger angle without discussing 
its implications, referring the reader to the thorough discussion 
presented in Ghisellini et al. (2002). 

The resulting line luminosity of the element X is then given 
by: ^ 

= ^^Ci'ion^ ~ 0.25^^LAfter (6) 
47r ?7cont 

where r]x is the luminosity ratio plotted in Fig. [l|. Note that 
the geometric factor C, is not relevant in this expression. For 
low ionization parameters ^ < lO-' one has ?7cont ~ 0.1 and 
r]x ^ 10"^, while for iron one obtains jype ~ 10"^ since its ef- 
ficiency is strongly decreased by Auger auto ionization (see 
also Ross, Fabian & Brandt 1996). Light elements can repro- 
cess into Ka line emission a sizable fraction of the continuum, 
yielding luminosities of the order of ~ 2.5 x 10'*'* erg s"' 
(Fx ~ 5 X 10"'"* erg cm"^ s"' at z ~ 1). In contrast, iron Ka 
emission will be one order of magnitude fainter, and therefore 
undetectable with present instrumentation. For larger ioniza- 
tion parameters 10-* < £, < 10"*, iron emission is efficient, while 
light elements are too ionized to emit line photons. A bright 
iron line should then be detectable in absence of light element 
features (see also Ballantyne & Ramirez-Ruiz 2001). 

3.2. Engine dominated models 

Evaluating the line luminosity in ED models is simpler and 
less affected by uncertainties. In fact, there is no geometric di- 
lution of the photons and the ionizing continuum is likely to be 
smaller than the afterglow luminosity at the time of the line ob- 
servation. Only the beaming factor f2/47r has to be taken into 
account, yielding a luminosity: 

n 

Lx = -r- ^After % ^ . 25 LAfter VX (7) 

47r 

In ED models, then, the line luminosity is a factor ?7cont less 
intense than in GD models. Lines from light elements (Mg, 
Si, S, Ar and Ca), should be very difficult to detect since 
i-IAr.Ca} < lO'*^ erg S-' {Fx ^ 10"'^ erg cm'^ s"* at z ~ 1). Iron 
lines should instead be detectable, in particular if the ionization 
is not too high and the spectral slope is flat (Ballantyne et al. 
2002). 

4. SUMMARY AND DISCUSSION 

The presence of emission features in the early X-ray after- 
glows of GRBs provides important clues for identifying the 
burst progenitor. We computed the luminosity of emission 
features from astrophysically relevant elements (Mg, Si, S, Ar 
and Ca) in reflection models and applied the results to the de- 
tectability of those features in the early afterglows of GRBs. We 
find that the line luminosity of these elements can be a strong 
indicator of the geometry of the reprocessing material. In par- 
ticular, light element lines are detectable only if the ionizing 
continuum is not directly observable. This geometrical set-up 
is naturally accounted for in GD models, in which the line pho- 
tons reach the observer after being reflected by material away 
from the line of sight. In ED models, the ionizing continuum 
and the line are observed simultaneously and for this reason the 
strength of light elements lines is greatly diluted, making them 
undetectable. 



Recently, a possible detection of Ka emission lines from 
light elements in the early afterglow of GRB 011211 has been 
claimed (R02). Reeves et al. claim that their spectrum is incon- 
sistent with a reflection model. We should however be aware 
that the reflection model used to fit the data does not include 
emission from S, Ar and Ca (Ross & Fabian 1983). In ad- 
dition, they fit a pure reflection model, without including the 
power-law external shock emission. This power-law compo- 
nent must however dominate the continuum, since the broad 
band lightcurve decreases as a power-law in time (R02). Let us 
now consider the Ka line of sulphur at 2.6 keV. It is detected at 
a flux level of (1 ±0.3) x lO"*'* erg cm'^ s"', when the [0.2-10] 
keV continuum has a flux of ~ 2 x lO"'-' erg cm"^ s"' . This flux 
ratio j]s = 0.05 ± 0.015 can be compared to the lower left panel 
of Fig. and Eq. 0. Even if the beaming correction approaches 
unity, the detected flux ratio exceeds the maximum allowed in 
an ED model, regardless of both the value of the spectral slope 
and ionization parameter. In GD models, however, the luminos- 
ity of light element features can be ten times larger due to the 
very small fraction of the continuum that is reflected for small 
ionization parameters (see the upper panel of Fig. H). 

Another important particularity of the observation of light 
element features in GRB 011211 is the lack of an iron line. 
If we compare this observation to that of the iron Ka line in 
GRB 991216 (POO), we find that the ratio of iron to sulphur lu- 
minosities is significantly different. The line ratio is Ls/^Fe < 
0.4 in GRB 991216 and Lg/Lpe > 2 in GRB 011211. This 
difference can be explained, as a consequence of the different 
ionization state of the reprocessing material - in GRB 011211 
the reprocessing material had an ionization parameter between 
30 < C < 300, while in GRB 991216 10^ < C < 10^ In the 
framework of GD models, the different ionization parameter 
can be interpreted as a different age of the SNR at the time 
of the burst explosion. If the remnant itself has a density that 
scales with the radius (and hence with age) hsnr oc ^snr' 
ionization parameter should also scale with radius for a constant 
ionizing luminosity. Young remnants should then produce light 
element lines, possibly in coincidence with weak nickel and/or 
cobalt lines, while old remnants should produce high ionization 
iron lines, and no light element features (see also Vietri et al. 
2001). In addition, the variability time scale of compact rem- 
nants should be smaller All these predictions are confirmed by 
observations - the features observed in GRB 01 121 1, in which 
a marginal evidence of Ni line is present, are observed to dis- 
appear after - 4 h (R02), while the iron Une in GRB 991216 
is detected ^ 18 h after the burst (POO). Is should also be re- 
marked that these conclusions do not depend on the assumption 
of a constant temperature of the reprocessing material. In fact, 
should this assumption be wrong or the temperature be esti- 
mated incorrectly, the result would be an horizontal shift of the 
lines in Fig. |l|. Since we have not discussed in which conditions 
a certain value of the ionization parameter is relevant, this does 
not affect the generality of our conclusions. 

In conclusion, the detection of Mg, Si, S, Ar and Ca lines 
with equivalent widths of hundreds of eV in the early after- 
glow of GRB 011211 can not be accounted for in any version 
of ED models. A two step explosion is necessary, in order to 
place the reprocessing material at a distance large enough to ex- 
plain the duration of the line and to let the ionizing continuum 
fade (hidden by the afterglow emission) before the lines can be 
observed. This evidence adds to the transient absorption fea- 
ture detected in GRB 990705 (Amati et al. 2000; Lazzati et al. 
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2001a) which strongly favors a GD model. Alternatively, the 
reflecting material may have been expelled during the common 
envelope phase of a binary progenitor system composed by a 
compact object and a helium star (Bottcher & Fryer 2001). In 
this case, all the discussion about GD models hold true, but a 
Ni line is not expected since the reprocessing material is not a 
SN explosion. 
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Fig. 1. — Ratio of the Ka line to the [1-10] keV continuum luminosities for various astrophysically relevant elements as a function 
of the ionization parameter ^ = 47rF[i_io]/n. Solid and dashed lines correspond to a power-law index a= 1.25 and a = 0.75 of 
the ionizing continuum, respectively. Black and grey lines correspond to solar and ten times solar metallicity of the reprocessing 
material, respectively. 
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Fig. 2. — Upper panel: ratio of the reflected to incident continuum in the [1-10] keV band for a solar metallicity slab. Two incident 
spectral slopes are considered, as in Fig. |l]. Lower panel: ratio of line luminosities of Mg, Si, S, Ar and Ca to the iron 
luminosity. Even though the iron line is generally more intense than all the other lines, there is a range of ionization parameters 
^ - 10^ for which iron Kq^ emission is quenched by Auger auto ionization, and light elements dominate. This figure is made for a 
continuum power-law index a = 1.0 and solar metallicity. Changing these parameters can be a moderate effect on the shape of the 
curves (see Fig. |l]). 



